Unification models with reheating via Primordial Black Holes 
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We study the possibility of reheating the universe through the evaporation of Primordial Black 
Holes created at the end of inflation. This is shown to allow for the unification of infiation with 
dark matter or dark energy, or both, under the dynamics of a single scalar field. We determine the 
necessary conditions to recover the standard Big Bang by the time of nucleosynthesis after reheating 
through black holes. 
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I. INTRODUCTION 

The quest for the unification of inflation, dark matter, 
and dark energy, in different combinations, by a single 
field, has been studied in Refs. p]-[l5|. The main mo- 
tivation behind all proposals is that we do not yet un- 
derstand the nature of the components responsible for 
the three phenomena, but we do know that their special 
properties are beyond the realm of the ordinary matter 
described by the standard model of particle physics. 

An extreme, most economical, possibility is that all 
three phenomena can be explained by the existence of 
one single field. As was first put forward in Ref. [g, the 
simplest option at hand is a scalar field (j) with a poten- 
tial of the form V{(t)) = Vq + (l/2)m^0^. The energy 
scale Vb is to be set at the tiny value of the observed cos- 
mological constant considered in the concordance ACDM 
model, and the mass scale of the field, m ~ 10~^mpi, is 
determined by the amplitude of primordial perturbations 
generated during inflation. 

A crucial ingredient in unification schemes that include 
the inflaton field is a proper and smooth transition from 
an inflationary era to the Hot Big Bang phase we arc all 
familiar with. There needs to be a reheating phase in 
which the infiaton field can transfer most of its energy 
to other relativistic particles, yet be able to survive the 
reheating process and attain the right amplitude to later 
become a subdominant dark matter/dark energy compo- 
nent in the radiation-dominated era. 

The viability of some phenomenological models of re- 
heating for unification models was explored in Refs. [2|,@1 ■ 
Some other, more detailed proposals have been explored 
in the literature, from extra periods of inflation at low 
energy scales 0], to the more popular gravitational par- 
ticle production at the end of infiation [e.g.[ll.[lll [l6l - [l8{ . 



This reheating mechanism is largely ineffective, however, 
because the relative proportion of radiation to scalar field 
Prci/ Pcfi evolves at the same rate as the proportion of grav- 
itational waves to scalar field pcw/p^ uM- ^^ ^^e case 
of steep potentials, inflation is followed by a kination era 
(a possibility first explored in the context of electroweak 
baryogenesis [20[). If the kinetic energy of the scalar field 
dominates for some time, the relative density of gravi- 
tational waves created during infiation can grow signifi- 
cantly and violate the observational bounds. Specifically, 
the process of Big Bang Nucleosynthesis (BBN) imposes 
a condition on the energy density of gravitational waves 
PGW; BBN requiring that [2l[ 



Pgw(^bbn) < 0.2pi.cl(iBBN) 



(1) 



In fact, the BBN bound in Eq. ([T]) has been used to con- 
strain and exclude some models of braneworld inflation 
and other uniflcation models j22| . 

In this paper, we explore an alternative mechanism to 
reheat the Universe, exploiting the evaporation of pri- 
mordial black holes (PBHs) produced at the end of the 
inflationary phase [23ll24| . and study its usefulness in uni- 
fication models. A small initial population of tiny black 
holes will grow significantly (in relative energy density) 
in a universe dominated by a stiff fluid. The evaporation 
of such black holes at the right time would produce the 
required amount of radiation more effectively than, for 
instance, the process of quantum particle production in 
an expanding space-time. 

In the PBH reheating scenario, the black hole compo- 
nent effectively reduces the relative density of gravita- 
tional waves as 



PGW 
PPBH 



^cnd 



(2) 



* |jhidalgo@aatro.una m.nix| 
t lurcna@fisica.ugto.rTDc| 
■f .a.liddleCisussex.ac.ukr 



Reheating takes place when the energy density ppbh is 
transformed into relativistic particles. Therefore, even if 
the densities pcw and ppbh were equivalent at the end 
of inflation, we would only require 3 e-folds of kinetic 



energy domination to meet the nucleosynthesis require- 
ments. PBH reheating therefore solves the problem de- 
scribed above. Another advantage of reheating through 
the evaporation of PBHs is that the unification field does 
not need to be coupled to other matter fields, hence it 
can later behave as a totally uncoupled dark matter com- 
ponent. PBH reheating for unificati on p urposes was first 
considered in Ref . ; see also Ref . [IJ] . 

Our aim here is to determine the general conditions 
under which the unification of inflation with either, or 
both, of the dark fields can achieve the appropriate re- 
heating with the aid of PBHs and set up the required 
conditions before BBN. 

The featured model of PBH reheating is described in 
detail in Sec. [Hi Together with this description, we estab- 
lish the conditions that the scalar field potential must ful- 
fil in order to undergo a successful PBH reheating. A par- 
ticular example of our proposal is illustrated in Sec. IHII 
Finally, in Sec. IIVI we frame our proposal in the context 
of unification models and draw some concluding remarks. 



II. PBH REHEATING 

There are three basic requirements for a successful 
PBH reheating: 

i) a significant number of PBHs produced at the end 
of inflation, 

ii) that PBHs become the dominant component in the 
Universe, and 

iii) the evaporation of the PBHs after they become the 
dominant component and before the time of BBN. 

The sequence described above introduces some tuning. 
We require the formation of black holes massive enough 
to last for a few e-folds while the Universe is dominated 
by the scalar field kinetic energy. Additionally, the same 
PBHs must evaporate before BBN, or before any prior 
well-tested process at higher energies in the standard Big 
Bang model. Each one of these conditions can be used to 
constrain the free parameters of the reheating process. 



A. Production of PBH and domination 

In recent years, several works have demonstrated the 
compatibility of the observed power spectrum of matter 
fluctuations at cosmological scales with the formation of a 
significant number of PBHs at the end of inflation (icnd), 
see e.g. Ref. [23I, or indeed due to preheating [23. More- 
over, the relative growth of ppbh with respect to other 
components is largest for PBHs formed at iend- This jus- 
tifles considering the black holes formed right after the 
end of inflation as those responsible for reheating. PBHs 
could also be originated by features in the inflationary 
potential. In the inflationary phase, these features break 



the slow-roll conditions momentarily, and may enhance 
the amplitude of fluctuations and black hole production 
at a specific scale after the end of inflation [e.g. [231. This 
presents an alternative to the production mechanism con- 
sidered here. 

To quantify the number of black holes at the end of in- 
flation, it is customary to deflne with /3 the initial mass 



fraction of PBHs at the time of formation tr 



The 



Press-Schechter formalism of structure formation pre- 
scribes that /3 is a function of the variance a of matter 
fluctuations. If we define P as the Gaussian probability 
distribution of d, with a its variance at that particular 
scale, then 



/3M=/m/ nS,a)dS, 



(3) 



where (5min — Wcnd l28|, and Wcnd = w^{Und) is the equa- 
tion of state at the end of infiation, which coincides with 
that of the inflaton field which is by then the dominant 
matter component. Here, /m ~ 3""^/^ is the fraction of 
the horizon mass which constitutes the black hole at the 
time of formation [23| . The amplitude of threshold inho- 
mogeneities Smin is determined by the equation of state of 
the dominant energy component at the end of inflation, 

WondQ 

The energy density of the produced PBHs evolves as 
pressureless matter (i.e. dust, with wpbh ~ 0), so that 
after time iend, the ratio of the PBHs to the scalar held 
energy densities is 



Ppbh 



/3(a) fait) 



1- (3{a) Vflcnd 



A(W) 



/3(^) 
l-/3(a)' 



=A(N) 



(4a) 



N 



A{N)^3 w^{N')dN', 
Jo 



(4b) 



where N is the number of e-folds elapsed from tend, cr 
is the mean amplitude of matter fluctuations, and A(A^) 
is the (exponential) growth factor of the PBHs. Eq. (|4]) 
takes into account the fact that typically the scalar held 
equation of state does not remain constant after the end 
of inflation. 

The mass fraction of PBHs will always increase as long 
as A(A^) > 0, that is: 



w^ > 0, 



for some <cnd < t < ikii 



(5) 



where ikin is the time when the kination era of the scalar 
fleld ends and it starts to behave like dark matter or dark 
energy. For definiteness, we take the equation of state at 
the end of the kination era to be Wkin '■= U70(tkin) ~ 0, 



^ Note that we evaluate tOond at times shortly after the end of 
infiation, where w can have a value between !«;„£ = —1/3 right 
at the end of inflation, and the maximum at the kinetic energy 
domination ui = 1. Also, the precise value of ^min is dependent 
of the shape of the initial configuration 130(1 . 



and the total growth factor is A(A^kin)- So far, we are 
assuming that, as in most cases, the dark field-hke be- 
haviour of the inflationary field appears at the onset of 
oscillations of the field around the minimum of the scalar 
potential. This need not necessarily be the case, but the 
final results are insensitive to the precise details of the 
dark matter/dark energy transition. 

It is also clear, from Eq. ^, that for PBH domination 
we require 



ln/3((7) > -A(7Vki, 



(6) 



Whether this condition is satisfied will depend upon the 
scalar field model and the exact evolution of its equa- 
tion of state in the post-infiationary era. One can fore- 
see, though, that some models may fail to satisfy condi- 
tion (O. 

Finally, we note that in the PBH reheating process we 
do not expect a significant PBH merging rate at the phase 
of domination. This is because the value of /3 is initially 
much smaller than 1, and the density of PBHs grows 
predominantly because of the rapidly diluting scalar field 
energy density [31 1. 



B. The inflationary energy scale 

The dominant radiation at the start of the Big Bang 
is generated by the evaporation products of PBHs. De- 
pending on the initial mass of the black holes, the time of 
evaporation, iov, can be shorter than, equal to, or larger 
than the time at which cj) starts behaving as dark mat- 
ter/dark energy, ikin (the relation between these times 
can be derived directly from the characteristics of the 
scalar field model under consideration). 

For the sake of clarity, let us consider the case in which 
^cv = ikin- The evaporation time of PBHs is determined 
by their mass, which is directly related to the energy scale 
at the end of inflation, p^nd — Ksnd, 



-1/2 



MpBH = Jm o'^PendH^A = 47r%/3/ 



Pcnd 



I HI I 4 

771 



^pl : 



Pl 



(7) 

where rripi is the Planck mass. The time of evaporation 
is given by [32, Isj] 



sec. (8) 



iev=6.3xlO-41<i>-l(MpBH)f^^^^ 



Here $(A/pbh) is a function of the directly-emitted 
species of particles and takes a value of order 10 for the 
PBH masses of interest [e.g. |3J]. Combining Ecj^s. ([7|) 
and ([5]), we find that 

/ 1/4 \ -6 

tev = 2.3 X 10"-^^$"^(A/pbh) ^^ sec . (9) 

\mpij 

This time must also be the beginning of the standard Hot 
Big Bang (HBB), at some point before BBN. This last 



requirement imposes a bound on the initial mass of the 
black holes and, consequently, on Pcnd- Indeed, in view 
of Eq. ([9]), the condition tov < teBN = 1 sec. implies 



K'nd > 1-6 X 10-"$i/^(MpBH)r 



ipi ■ 



(10) 



C. The right amount of radiation 

On the other hand, we can determine whether radi- 
ation domination is achieved after the evaporation of 
PBHs. Assuming again that the standard HBB is recov- 
ered at the time of evaporation, we can write a simple 
expression for the ratio of PBH to scalar field densities 
as 

1/2 



/ PPBH \ 

V P4> A 



Prcl 
Pdm 



3 X lO 



26 



v: 



1/4 \ 3 



JTlpl 



fll) 



where we have considered that from tev onwards the Uni- 
verse is radiation dominated, and gev (ffeq) are the rcla- 
tivistic energy degrees of freedom at that time of evapo- 
ration (matter-radiation equality), respectively. 

Radiation domination is ensured as long as the evolu- 
tion of the scalar field equation of state allows the match- 
ing between Eqs. Q and pT|) at the time of PBH evap- 
oration. We will see that this is possible even when the 
constraints on the products of PBH evaporation are con- 
sidered. 



D. Additional constraints 

There are two important bounds on the density of the 
products of evaporation of PBHs. These can be written 
in terms of the energy density at the time of black hole 
formation and, therefore, introduce bounds to the unifi- 
cation inflationary potential. The first constraint comes 
from the consideration that a PBH could leave behind 
a Planck mass relic after evaporation. Such relics must 
not exceed in density the dark matter component [2ll . l34| . 
i.e., we require 



/^rclic 

Pdm 



<1 



(12) 



From Eq. ([7|) we can read the proportion between the 
PBH mass and the Planck mass of its relic. Therefore 
the energy density ratio of PBHs with relics is 



/ PPBH \ _ / Prcl 
\ Prclic / cv V Prclic 



47rV3/; ' ' ™P' 



TM 



V„ 



1/4 



(13) 



Combining Eqs. pTJ) and ([T3|), the condition ([T2|) is met 
when 



V}J* < 6.8 X 10- 



end 



V- 



(14) 



If the condition p2)) is met after evaporation, then the 
same ratio of rehcs to dark matter is maintained at sub- 
sequent times. 

The combination of Eqs. (|T0|) and ([M]) leaves but a nar- 
row window for the range of energies at which inflation 
must end if we want the PBHs to reheat the Universe. 
The tight constraint can however be relaxed if no Planck 
mass relics are left over after evaporation [35|. 

This is a crucial point for standard (slow-roll) inflation- 
ary scenarios, for which the end of inflation takes place 



at energies V^^^ 



10 



-3, 



ipi. PBH reheating in this case 



is only possible if no Planck relics are left over after evap- 
oration, or if the black holes responsible for reheating are 
created at energies below p^nd- 

A second constraint comes into play in the form of 
gravitational waves. When black holes evaporate, a con- 
siderable fraction of energy is emitted in the form of 
gravitons. Paradoxically for the initial motivation of this 
paper, it appears that in the case of a dominant energy 
density component of PBHs, a significant amount of grav- 
itational waves is produced at evaporation. This could 
violate the very constraint we intend to alleviate. 

When PBHs dominate and then evaporate, the amount 
of gravitational radiation produced is |3l| , 



r^GW = 0.36/.go 



(15) 



We are interested, however, in times of evaporation prior 
to nucleosynthesis, when ^cv > 10.75. Consequently, 



V Prcl 



f^Gw(tcv)<0.03, 



(16) 



which readily satisfies the gravitational waves bound in 
Eq. HD. 



1 X 10" 




1 X 10" 



1 X 10" 



1 X 10" 



1 X 10" 



III. EXAMPLES 



FIG. 1. The curves indicate tlie required variance of the 
matter spectrum a^, see Eq. Q, to create enough PBHs 
that subsequently evaporate and account for all the radia- 
tion in the early Universe before BBN. The figure at the 
top considers pcnd = (10~ rripi) and the bottom figure takes 
Pcnd = (10~'^mpi)*. We plot the scales of the evaporation of 
PBHs and of BBN with a grey and a black line, respectively. 
The area below pbbn is excluded by observations. 



To illustrate our proposal, let us look at some scenar- 
ios where PBH reheating can take place. In the model 
considered in Ref. |2|, which intended a unification into 
a single scalar field of inflation and dark matter only, as 
in many other braneworld models, inflation ends at an 
energy scale of V^^^ « 10~^TOpi, which is within the al- 
lowed values of constraints (|10p and ([M]) . The expansion 
of the Universe between the end of inflation and the onset 
of oscillations is such that 



flkir 
^cnd 



= 60; 



(17) 



and the scalar field enters a kination phase in between, 
during which w^ ~ 1 . This means that the growth factor 
in Eq. (g]) is A(Arkin) ^ 12. In consequence, Eq. (O teUs 
us that to reach the domination of PBHs before (j) starts 
oscillating, we must havc/3 > 10^^. This is in agreement 
with the results in Ref. [2|. 



Another example is the quintessential infiation model 
in [l| , for the unification of infiation and dark energy. Af- 
ter inflation, the field goes into a kination era for which 
the Universe expands akin/oond = 10* ~ e^*. Thus, the 
quintessential model can be more easily implemented, as 
it only requires /3 > 10~^^ to generate enough radia- 
tion from PBHs. Moreover, if the inflation ends due 
to an instability of the hybrid infiation kind, then the 
smaller-scale fiuctuations would show an enhanced am- 
plitude [S^, rendering a larger probability of PBH for- 
mation. 

In Fig. [T] we present graphically the valid parameter 
values for a successful PBH reheating. The curves in 
both figures indicate the mean amplitude of matter fluc- 
tuations, a in Eq. (|3]), required for a given energy density 
(or time) at which the scalar field starts oscillating. 

The figures show the importance of the difference be- 
tween the equation of state right after the end of in- 



flation Wend and at the subsequent kination epoch w^. 
For simphcity, we have assumed that w^ = const, dur- 
ing the kination phase of the scalar field. This is not a 
strong restriction, as the total growth factor can be writ- 
ten in terms of a mean equation of state as A(A'kin) = 
3wmoan-/Vkin, and tlicn the plots in Fig. [1] can be inter- 
preted in terms of lUmoan too. 

In particular, we see that the small value Wond < 0.1 
allows the formation of large populations of PBHs with 
a relatively small variance. Note that the potential of 
the figure at the bottom violates the constraint imposed 
by Planck relics, a constraint that, as mentioned above, 
need not necessarily apply. 

By considering a range of values for pkin above (below) 
the evaporation scale pov in the top (bottom) figure, we 
are showing that the conditions required for PBH reheat- 
ing when tov 7^ ikin are not too different from those in 
the case of equality (described in the previous section) . 

For the sake of generality we are not computing here 
the amplitude of fluctuations for specific models at the 
end of inflation. It suffices to mention that for some 
models with Wcnd ~ 0, and under suitable conditions, the 
production of a considerable amount of PBHs is possible 
at horizon and subhorizon scales [33, S3 ■ 



IV. UNIFICATION MODELS AND FINAL 
REMARKS 

If the unification of infiation with dark matter or dark 
energy, in any given combination, is to be achieved by a 
single field, then an efficient reheating process is needed 
which docs not heavily rely upon the decay of the unifi- 
cation field. 

We have seen that, in the presence of a large enough 
population of PBHs, unification models can be possi- 
ble if the unification field enters a long enough kinetic- 
dominated era, so that the PBHs come to dominate the 
matter budget. PBH domination is easier for long kina- 
tion periods and stiffer values of the effective equation of 
state w^, but the latter also requires a larger value of the 
mean amplitude of primordial perturbations a at the end 
of infiation. 

In this article we have not included an explicit mecha- 
nism for the PBH formation, as would be required for a 
fully self-contained model. Most economical, if possible, 
would be for the unification field itself to feature per- 
turbations rapidly growing at late times, as may happen 
if its second derivative becomes large as infiation ends. 
Alternatively, a phase transition (for instance of hybrid 
inflation form) may lead to enhanced perturbations as it 



takes place [e.g.[2J, [3a,|33. In any case, it is worth men- 
tioning that the mechanism required for PBH formation 
need not spoil or interfere with the dynamics followed by 
a unification model in the dark matter or dark energy 
eras. 

Because PBHs form right after the end of infiation, 
most of the observational restrictions can be written in 
terms of the energy scale at the end of infiation. In the 
simplest PBH scenario, the constraints are satisfied only 
if this energy scale is a few orders of magnitude below the 
Planck value. One should bear in mind that the possibil- 
ity of production of Planck relics renders PBH reheating 
incompatible with standard inflationary models, unless 
the energy scale at the end of inflation happens to be of 
the order of 10~^ inp\. 

The strongest constraint on the featured mechanism 
comes from the evolution of the scalar field equation 
of state after inflation. The unification model should 
have a long enough kination period to allow PBH 
domination before evaporation. Typical examples are 
braneworld inflation with an exponential potential [2, |40[ 
and quintessential inflation [l|. However, even if the 
Planck relic constraint could be set aside, this is a con- 
dition that many inflationary models will not be able to 
surpass. 

Our results rely heavily on the assumption that re- 
heating proceeds from PBHs formed right at the end of 
inflation. The results can change if PBHs are formed af- 
ter the end of inflation, or if we consider the kination 
period to happen beyond BBN (a possibility explored in 
Ref. [41|). Therefore, the parameters of the theory can 
be flexible beyond the constraints presented here. 



ACKNOWLEDGMENTS 

J.C.H. gratefully acknowledges sponsorship from 
DGAPA-UNAM and the support of PAPHT-UNAM 
(grant IN116210-3). A.R.L. was supported by the Sci- 
ence and Technology Facilities Council [grant number 
ST/1000976/1]. L.A.U.-L. thanks the Berkeley Center 
for Cosmological Physics (BCCP) for its kind hospital- 
ity, and the joint support of the Academia Mexicana 
dc Ciencias and the United States-Mexico Foundation 
for Science for a summer research stay at BCCP. This 
work was partially supported by PROMEP, DAIP-UG, 
and by CONACyT Mexico under grants 56946, 167335 
and 10101/131/07 C-234/07 of the Instituto Avanzado de 
Cosmologia (lAC) collaboration. 



[1] P. J. E. Peebles and A. Vilenkin,'P hys.Rev.| D59, 063505 

(1999), arXiv:astro-ph/9810509 
[2] J. E. Lidsey, T. Matos, and L. A. Urena-Lopez, 
,Phys. Rev., D66, 023514 (2002), ,arXiv:astro-ph/0111292, 



[3] T. Padmanabhan and T. R. Choudhury,' Phys.Rev.| D66, 

081301 (2002), arXiv: hep-tli/020505 5 
[4] R. J. Scherrcr, Phys^RevXetEf 93, 011301 (2004), 

iarXiv:astro-ph/0402316] 



[5] 

[6] 
[7] 

[8] 



[9] 
[10] 

[ii: 

[12] 

[13] 
[14] 

[15] 

[16] 
[17] 

[18] 



[19] 

[20] 
[21 



T. Matos, J.-R. Luevano, and H. Garcia- 
Compean, |Int.J.Mod.Phys. A23, 1949 (2008), 



[22] 



[23] 



103512 
(2007), 
127301 



|arXiv:hep-th/0511098| 

A. Arbey, IPhys-RevT] D74, 043516 (2006), 

|arXiv:astro-ph/0601274r "^ 

V. H. Cardenas, Phys.Rev. D73, 
(2006), arXiv:gr-qc/0603013 D75, 083512 
|arXiv:astro-ph/0701624G. Panotopoulos, D75, 
(2007), arXiv:0706. 2237 [hep-ph], 

A. R. Liddle and L. A. Urena-Lopez, |Phys. Rev. Lett.| 
97, 161301 (2006), 'arXiv:astro-pti/0605205A. R. Lid- 
dle, C. Pahud, and L. A. Urona-Lopez, Phys. Rev. D77, 
121301 (2008), arXiv:0804.0869 [astro-ph] 

C.-M. Lin(2009), arXiv:0906.5021 [hep-ph], 

A. B. Henriques, R. Potting, and P. M. Sa, |Phys.Rev.| 
D79, 103522 (2009), [arXrv:0903.2014 [astro- ph.COp 
N. Bose and A. S. Ma,junidar, |Phys.Rev.| D80, 103508 



(2009), |arXiv:0907.2330 [astro-ph. CO]' 
A. de la Macorra and F. 
arXiv:0910.163 2 [astro-ph.CO 
arXiv:0911.5172 [astro-ph.CO 



Briscese(2009), 



;2009), 

J 



M. Bastero-Gil, A. Berera, B. M. Jackson, and A. Taylor, 
|Phys Lett.| B678, 157 (2009),'arXiv:0905.2937 [hep-ph]J 
R. N. Lerner an d J. McDonald, Phys.Rev. D80, 123507 
(2009), arXiv:0909.0520 [hep-ph] M. Okada and Q. Shafi, 

D84, 043533 (2011), arXiv:1007.1672 [hep-ph], 

J. De-Santiago and J. L. Cervantes-Cota, |Phys.Rev.| 
D83, 063502 (2011), [arXiv:1102.1777 [ astro-ph.COp 
L. H. Ford and L. Parker, |Phys. Rev.| D16, 1601 (1977) 

B. Spokoiny, 'Phys. Lett.' B315, 40 (1993), 

|arXiv:gr-qc/9306008 [gr- qc] 

G. Hue y and J. E.n:. idsey, P hys.Lett.j B514, 217 
(2001), |arXiv:astro-ph/0104006^. Dimopoulos and 
J. W. F. VaUe, [Astropart. P^ 18, 287 (2002), 



|arXiv:a stro-ph /0111417[ 

V. Sahni, M. S ami, and T. Souradee p, |Phys. Rev\] D65, 

023518 (2001), arXiv:gr-qc/0105121 

M. Joy ce, Phys.Rev^ ^D55, 1875 (1997), 



|arXiv:hep- ph/9606223 [h ep-ph] 

V. F. Shvartsman, Pisma Zh. Eksp. Teor. Fiz. 9, 

315 (1969)B. J. Carr, Astronomy and Astrophysics 

89, 6 (Sep. 1980)R. H. C yburt, B. D. Fiel ds, K. A. 

Olive, and E. Skillma n, [Astropart. Phys.| 23, 313 

(2005), arXiv:astro-ph/0408033B. J. Carr , K. Kohri, 

Y. Sen douda, and J. Yokoyama,,Phys.Rev.| D81, 104019 

(2010), |arXiv:0912.5297 [astrcHph.CO]| 

M. Gi ovannini, |Phys. Rey^ D58, 083504 

(1998), |arXiv:hep-ph/9806329p A. Boyle and 

A. Buonanno, Phys.Rev. D78, 043531 (2008), 

arXiv:0708.2279 [ast ro-ph] 

J. D. Barrow, E. J. Copeland, E. W. Kolb, and A. R. 

Liddle, |Ph^s^ev]D43, 984 (1991) 



[24] 
[25] 



[26] 
[27] 
[28] 

[29] 
[30] 

[31] 

[32] 
[33] 

[34] 

[35] 

[36] 
[37] 



[39] 

[40] 
[41] 



J. Garcia-Bellido, A. D. Linde, and D. Wands, fPhys.RevT] 

D54, 6040 (1996),,arXiv:astro-ph/9605094 

IJCAPI 0701, 
Kohri, 
, 038 
and 
(2009), 
Sasaki, 



S. Chongchitnan 
Oil (2007), 
D. H. Lyth, 
(2008), _ 
K. Kohri, 



and G. Efstathiou, 
arXiv:astro- ph/0611818k. 
and X! 



Melchiorri, 0804 

arXiv: 0711.5006 [hep-ph] L. Alabidi 
D80, 063511 
arXiv:0906.1398 [astro-ph. CO]L Zaballa and M 
JCAPI 1003, 002 (2010), | arXiv:0911. 2069 [ast ro-ph.CO]] 



A. M. Green and K. A. Malik |Phys. Rev.[ D64, 021301 
(2001), |arXiv:hep -ph/0008113' 

P. Ivanov, Phys.Rev. D57, 7145 (1998), 
|arXiv^ stro-ph/9708224 [astro-ph] 

Carr, Astrophys. J.. 201, 1 (1975)1. Musco, J. C. 
and L. RezzoUa, [Class. Quant. Grav.l 22, 1405 



B. J. 

Miller, 



L. RezzoUa, [Class. Quant. Grav.| 22 

(2005), [arXiv:gr-qc/0412063' 

B. Carr, J. Edited by Sanz, and L. Giocoechea, m Ob- 
servational and Theoretical Aspects of Relativistic As- 
trophysics and Cosmology. Proceedings, International 
Course, Santander, Spain, September 3-7, 1984 (1985) 
A. M. Green, A. R . Liddle, K. A. Ma- 
lik, and M. Sasaki, [Phys.Rev.[ D70, 041502 
(2004), larXiv: astro-ph/0403 1 8 IX C. Hidalgo 

and A. G. Polnarev , D79, 044006 (2009), 

arXiv:0806.2752 [astro-ph] 

R. Anantua, R. Easth er, and J. T. Giblin,iPhys. Rev. Lett. [ 
103, 111303 (2009), arXiv:0812.0825 [astro-p"h]] 

S. Hawking, Co mmun.Math.Phys.[ 43, 199 (1975) 

H. L Kim, C. H. Lee , and J. H. Ma cGibbon, [Phys.Rev!] 

D59, 063004 (1999), ■arXiv:astro-ph/9901030 

A. S. Josan, A. M. Green, and K. A. Malik, iPhys.Rev.[ 
D79, 103520 (2009), arXi v:0903.318 4 [astro- ph!CO]] 
S. D. Hsu, ,Phys.LetE] B644, 67 ^2007), 
|arXiv:hep-th/0608175p ! Carr, L. Modesto, and 
I. Premont-Schwarz(2011), |arXiv:1107.0708 [gr-qcj] 
D. Roberts, A. R. Liddle, and D. H. Lyth, [Phys.Rev^ 
D51, 4122 (1995), arXiv:astro-ph/9411104 [astro- ph]|' 
D. Lyth, K. A. Mal ik, M. Sasaki, and I. Zab alla, JCAP [ 
0601, Oil (2006), |arXiv:astro-ph/0510647| r. Suyama, 
T. Tanaka, B. Bassett, and H. Kudoh, 0604, 
001 (2006), [ arXiv:hep-ph/0601108f <. Jedamzik, 
M. Lemoinc, and J. Martin, 1009, 034 
arXiv: 1002.3039 [astr o-ph.CO] 
JCAP 



1107, 



035 



(2010), 
(2011), 



, 028 
and 



(2011), 

P. Kli- 

H. 



D. H. Lyth, 
|arXiv: 1012.4617 [astro-ph.CO], 

E. Bugaev and P. Khmai, IJCAPl 1111 

'arXiv: 1 107.3754 [astro -ph. CO] E. Bugaev 

mai(2011), [arXiv :1112.5601 [astro-ph.CO], 

Lyth(2012), arXiv: 1201. 4312 [astro-ph.COj"^ 

E. J. Copeland, A. R. Liddle, and J. E. Lidsey, [Phys.Rev.] 
D64, 023509 (2001), arXi v:astrcHph /0006421 

S. Dut ta and R. J. Scherrer, [Phys.Rev., D82, 083501 
(2010), arXiv:1006.4166 [astrc^ph.CO]] 



